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We derive an expression which can be used in a relatively straight- 
forward manner to obtain either the autocorrelation function or the spectral 
density of any "reasonable" function of any stationary gaussian process. 
The expression is used to study the spectral density of a sinusoidal wave 
which is phase modulated by a hard-limited gaussian process. 

I. INTRODUCTION 

A band-limited gaussian process having a rectangular power spec- 
trum is assumed for some purposes to be an adequate approximation 
for several classes of modulating signals encountered in phase modu- 
lation (PM) systems. 1 However, in an actual implementation of a PM 
system, the modulating signal passes through circuitry which saturates 
when the signal rises above a fixed level. This clipping (i.e., limiting) 
level is usually adjusted fairly high (nominally at four times the rms 
of the modulating signal) and then ignored in any subsequent analysis 
of the system. Consequently, the objective of this study is to deter- 
mine the qualitative effect of hard-limiting the modulating signal in 
a PM system. 

From a mathematical viewpoint, we can obtain an understanding 
of the preceding question by investigating the following problem: 
Find the spectral density of a sinusoidal wave which is phase modu- 
lated by a function g(X t ) of the stationary gaussian process X t . Of 
course, this version of the problem can also be viewed as finding the 
spectral density of a (composite) nonlinear function of a gaussian 
process; a problem originally studied by S. 0. Rice, 2 D. Middleton 3 
and W. R. Bennett. 4 In fact, we do use their approach (representing 
the nonlinearity in terms of a transform) to derive an expression 
which is essentially the starting point of our analysis. However, using 
our relation avoids some of the complexity associated with the trans- 

1025 



1026 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1971 

form method. A derivation of the expression is given in Appendix A. 
The notation and general results are presented in Section II along 
with two examples. In the third section, we obtain specific results for 
the hard-limiting case. 

II. SPECTRAL DENSITY OF A PM WAVE MODULATED BY A NONLINEAR 
FUNCTION OF A GAUSSIAN PROCESS 

Let W(t) denote a constant-amplitude sinusoidal wave which is 
phase modulated by a real- valued function g(X t ) of the stationary 
gaussian process X t . That is, 

W(t) = A cos ( Wc * + g(X t ) + B) (1) 

where A is the wave amplitude, f e = w c /2tt is the carrier frequency, 
and 6 is a random variable with probability density function 



T 9 (d) = 



£- for ^ 6 < 2ir, 

2X (2) 



.0, otherwise. 

To obtain the spectral density for W(t), it is convenient to express 
the wave in terms of complex variables as 5 

W(t) = Re {W(t)\ (3) 

where* 

W(t) = A exp \iiftjb + $)]V(t), (4) 

and 

V(t) = exp [jg{X t )\. (5) 

Since X t is stationary and 8 is uniformly distributed, both W(t) and 
V(t) are wide-sense stationary. 5 Moreover, the spectral density S w (j) 
of W(t) is given by 

SM = t[S.(j - U) + S v (-f - /.)] (6) 

where S v (f) is the spectral densi ty of V(t). 5 

If we define R v (t) = (V(t + r)V{t)), the Wiener-Khintchine theorem 
implies that 

S v (j) = f M R,(t) exp [-J27r/r] dr. (7) 

J — oo 



* We use Re (z) and Im (z) to denote respectively the real and imaginary parts 
of the complex variable z. The symbol z denotes the complex conjugate of z. 
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So, assume that X t + T and X t are gaussian with joint probability 
density function 

/ x 1 { x\ — 2rx^x 2 + x\ 

p(x 1 , x 2 ) = ———=== exp ' 



2 1 rTWl-r 2 "V 2r*(l-r a ) 

- oo < x,- < co , (8) 

where B x (r) = (X t+T X,), T 2 = R x {0) and r = fl, r ( T )/r 2 . It follows 
that 

/ exp [jigfa) - g{x 2 ))]-p{x l , x 2 ) dx 1 dx 2 . (9) 

• 00 w —oo 

In Appendix A, we show that if G(x) is of exponential order and 
of bounded variation on bounded intervals/ then 



/CO /iCO 00 

/ Gix^Gix^pix, , x 2 ) dx, dx 2 = 2 a n r n , \ r \ £ 1, (10) 
■oo "-co n = 

r 2n \ r ( d n \ 
a " = ^T \L G{x) \d? p{x) ) dx 



where 



and 



*-;7b-(-£) 



V27rr 
is the probability density for X,. Setting (7 (a:) = exp[;g(.r)], 

fl.w = E cy\ (ii) 

n-0 

where* 

p2n I »co T ,n II 2 

C " = ^7 I /.„ GXP ^(x)]\J^p(x)j dx I • (12) 

Recall that the Hermite polynomial of degree n, H„(x), satisfies the 
relation 

£- n exp (-x 2 ) = (-l) n H n (x) exp (-x 2 ). 



t We define such functions to be "reasonable." 

* Since |exp [jg(x)1\ = 1 for all finite g(x). it follows that equations (11) and 
(12) are valid for all functions g(x) which are of bounded variation on finite 
intervals. This should include most functions which one might encounter in 
engineering applications. 
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Consequently, 

Cn = 2^! | £, exp ^(x)W.(-^y)p(x) dx '• (13) 

Since the coefficients {C„} are independent of B a (r), an infinite 
series representation for the spectral density S v {f) can be obtained 
from equations (7) and (11). To obtain the series, let 



St(f) = f R z (t) exp (-jfcr/r) di 

•/ — to 



(14) 



and define S n (j) to be the n-fold convolution of 8i(f) with itself; i.e., 
S n (f) = 5,_i (/) * Si(J) for n ^ 2. It follows that 



S,(1) = C 8(f) +IC„^?, 

n = 1 1 



— oo < / < » . 



(15) 



Using equations (13) and (15), one can see that the carrier power 
is given by 



C = |<exp[;<KX ( )]>| 2 . 



(16) 



2.1 Examples 

To obtain a better understanding of the preceding results, we present 
two examples. The first example uses g(x) — x (and serves as a partial 
check on our results since this case is well known 5 ). The second 
example assumes extreme clipping, 

, . ! b for x > 0, 
gix) = \ 

{-b for x < 0, 

and serves as a preview for the hard-limiting case presented in 
Section III. 

2.1.1 Phase Modulation 

For this example, g{x) = x for all x. Consequently, from equation (13) 
we have, for n ^ 0, 



C« = 



1 
2 n n\ 



Thus, 






^2n 


1 


2 


2 2n (2n)! 


V, 



1 f °° 

-y= I exp (;V§ Tx)H n (x) exp (— x 2 ) dx 
fir *-« 

7= / cos (V§ rx)# 2n (x) exp (-x 2 ) dx 

7T •'O 



(17) 



(18) 
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C 2n = exp(-r 2 ) 



flpy 
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(19) 



(2n)! ' 

(see Ref. 6, Section 7.388-3) . 

In a similar fashion C 2 „ + i = exp(-r 2 ) (r 2 ) 2n+1 /(2n + 1)!, which 
implies that 



C„ = exp(-r 2 ) :L T for n ^ 0. 



m 



(20) 



Now, substitute equation (20) into equation (11) to obtain the well- 
known result 5 

R,(r) = exp[-r 2 (l -r)]. (21) 

2.1.2 Phase Modulation by an Extreme-Limited Gaussian Process 
For this case, 



M-\ " " x>0 ' 

[-b if x < 0. 
Using equation (22) in equation (13), one obtains 



(22) 



C --*- 
° n 2 B n! 



exp (jb) J HJ^^jpix) dx 

+ exp (-jb) J_ ff^^l^jpta;) da; 



(23) 



which reduces to 



c - = iwfef ff '" Wexp( - x!) 4 (24) 



and 



ft* - ^yfm (^ jf ^-W exp (-x') &)', (25) 
for n ^ 0. Since 

J" exp {-y 2 )H n {y) dy = #^(0) for n £ 1, 



we have 



Co = (cos (b)) 2 , 
C 2n = for n ^ 1, 



(26) 

(27) 
(28) 
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and 

**- s^Sflfc^W- (29) 

Using the fact that H (0) = 1 and # 2n (0) = (-l) n 2 n [l-3 • • • (2n - 1)] 
for n ^ 1, we see that 

C x - - (sin (6)) 2 

7T 

and 

C -« " i ** <»>' VZ Tn^n 2 Zl) l) fOT " * L (30) 

Substituting equations (27), (28) and (30) into equation (11) and 
recalling the series expansion for Arcsin x, one obtains 

R,(r) = (cos (6)) 2 + (sin (b)) 2 \ Arcsin (|^) , (31) 

for — 0° < T < oo . 

Notice that the carrier power, Co = (cos (6)) 2 , vanishes if b = ir/2-\-kir, 
k = 0, 1, • • • , while all the power goes into the carrier whenever b = kir, 
k — 0, 1, • • • • If b J* kir, the continuous part of the spectrum is simply a 
scaled version of the spectrum for an amplitude-modulated wave when 
the modulator is an extreme-limited gaussian process. This problem was 
studied by J. H. VanVleck. 7 

III. PHASE MODULATION BY A HARD-LIMITED GAUSSIAN PROCESS 

A problem of interest arises when the function g represents an ideal 
hard-limiter; i.e., 



g(x) = 



b for x > b, 

x for | x | ^ 6, (32) 

— b for x < — b, 



for some b ^ 0. 



3.1 Carrier Power 

Noting equations (13) and (16), one can see that the carrier power 
is given by 
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x- ' ' 2 



Co = ( vtp L cos {x) exp (~w) dx 

Defining a = b/T (the relative clipping level), we have 

C = (-^= j ' cos (\/2 Tx) exp (-x 2 ) dx + cos ib) erfc MgJJ (33) 
where 

erf (£) = —= [ exp (-t 2 ) dt (34) 

is the error function 8 and erfc (z) = 1 — erf (2). 

A more useful expression for C» can be obtained in terms of the 
error function with complex argument. Apparently, 



9 «a/vz 

—7= / cos (\/§ Tz) exp (-x 2 ) dx 

Vir •'0 

= Re S-V [ exp (-i\/§ Tx- - x 2 ) dx\ , 

= exp (-|) Re {^= /_"' ^ exp (-(, + j ±J) dx 

= 1 exp (-£) Re {erf (^ + i£)- erf fe + i£)}, (37) 



(35) 
(36) 



r\ _ j ,/ a , . r v 



= exp^- T jRejerf^ + ;^. ;; , 

Consequently, 

Co = (exp (-£) Re {erf (^ + j ^)} + cos (b) erfc (^))' (39) 

can be calculated numerically. s However, we can obtain a better 
analytical understanding of the carrier power by expressing the 
integral in equation (36 > as 



+ To make the last step, we used the relations 8 



erf (-2) = -erf (z) and erf (2) = erf (2). 
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-p (-t){ rc ^ r" exp (-( x + 4)1 4 

Since exp (— z 2 ) is an entire function, the integral in equation (40) is 
independent of the path taken from i(r/V2) to a/v2 + i(r/V2~). It 
is useful to use a contour consisting of three straight line-segments: the 
first from j(r/vl) to (on the imaginary axis), the second from to 
a/VZ (on the real axis) and the last from a/V2 to a/v2 + j(T/y/2). 
Integrating along these lines, it is straightforward to show that 

expl— tt) ReS— 7= / , exp (-z 2 ) dz 

\ II LV7T JHT/V2) 

= exp (-£) erf (^ 

+ exp (-£) -^= jT'' "am (vS ay) -exp (-(y - ?/)) 4y (41) 
Combining equations (33) and (41) yields 

Co = (exp (-£) erf (^) + cos (b) erfc (^) 

+ exp (-|) ^ f^sm^ay) exp(-(y - /)) %)'• (42) 

When r « 1 (low-index case), the last term inside the brackets is 
negligible so that 

C tt (exp (-y) erf (^) + cos (b) erfc (^J when T « 1. (43) 

It is interesting to observe that erf (a/V2) is the probability that X, 
is inside the clipping levels [ — 6, b] [and erfc (a/V2) is the probability 
that X, is outside the clipping levels]. Hence, we can reason that 
exp (— r 2 /2) is the average dc voltage associated with V(t) when X t is 
inside the clipping levels while cos (6) is the average voltage when X, 
is outside the limits. [One can see from equation (16) that the carrier 
power Co = | (VU)) | 2 .] 

Pursuing these thoughts further, one might conclude that the 
integral in equation (42) represents a high-order correction (to the 
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preceding approximation) associated with transitions of X t across the 
clipping thresholds. When T is small, these transitions are rapid and 
the correction (i.e., the integral) is small. However, as T increases, the 
transitions are slower (at least when X, is bandlimited) and the 
integral in equation (42) can increase in magnitude for certain values 
of a. Moreover, 

— j= j sin (V2" ay) exp (-(— - VJ) dy 

« sin (J>)F\^j when r » V2~, 



where F(x) = 2/Vtt /S exp [-(.r 2 - t 2 )] dt is Dawson's function. 8 
Consequently it follows that when r » v2 (i.e., the high-index case), 

Co « (exp (-£) erf (^) + cos (6) erfc (^) 

+ sin(6)exp(-|)F(^)) 2 - (44) 
Thus, when a is not too large, 

Co « (cos (b) erfc (^) + sin (b) exp (-t)f(^JJ ■ (45) 

The quantity inside the brackets is a sinusoidal function of a. Hence, 
there are various values of a (for fixed D which will completely null 
the carrier power C„ while other clipping levels give rise to a relative 
maximum carrier power. 
Using 

(V(t)) = exp (-|) Re {erf (^ + j ^)} + cos (b) erfc (^) , 

we computed exp (T s /2)(V(t)) as a function of a for several values of r. 
Some of the results are displayed in Fig. 1. (Since the dynamic range of 
this function is so large, we used a linear scale between — 1 and 1 on the 
ordinate axis and a logarithmic scale otherwise.) 

The general features of the three functions are as indicated above. 
Since (7(0) = 1 whenever a = 0, each curve starts at exp (r 2 /2) 
for a = 0. For T ^ 1 (not shown in Fig. 1) the curve decreased mono- 
tonically to one. For T = 2, this relative average remains positive but 
oscillates slightly before converging to unity. However, for the high- 
index case, the relative mean is similar to a damped sinusoid (as dis- 
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Fig. 1 — Relative mean of V(t) as a function of relative clipping level, a is 
the ratio of the clipping level b to the rms modulation-index r. 



cussed above). This extreme oscillation about zero apparently occurs 
for T ^ 3. Notice that the value of a required for the relative mean 
(and therefore the carrier power) to stabilize increases as T increases. 
Hence, for high-index modulation, setting a W 4 does not necessarily 
guarantee that the carrier power is unaffected by the clipping. 

3.2 Continuous Part of the Spectrum 

Numerical techniques are required to obtain the continuous part of 
the spectrum. However, one can see from equations (23) and (31) 
that the coefficients {C„} in the series are independent of R x (t) and 
need only be computed once for any particular configuration of the 
hard limiter (and choice of r). Consequently, it seems worthwhile to 
derive an algorithm for efficient computation of {C„}. 

Noting equations (11), (12) and (13), one can see that 



R.(r) = S 



R(t) 



c n 



(46) 



where 



C„ = 



r 2n I !»(») 



ni 



(47) 
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and 

JXn) = /" exp [-jg(x)](£i P(x)) dx. (48) 

Consequently, 

- / exp (-jx) dp(x) + exp (-jb) / dp(x) + exp(;6) / dp(z) 

J-b J b •'-co 

= exp(-|)Re(erfc(^ + 4)). (49) 

Next, for n 2: 2 we have 
j-V(n) = / exp (-jr) ^^np(a;)j + exp (-#) J d\^=rp(x)J 

+ exp(;6)/_ oo d(£^ P (x))- 
Integrating the above expression by parts, we have 

v(n) = y=t J exp (-j'z)(^prrp(x)J dx. (50) 

Recall that 

£- **» (-If) = iwr) H iwf) exp \-w) 

so that integrating equation (50) by parts obtains 

V (n) - * - 1) - ^j (^) Hn-a(^r) exp (-jr - ^) [^ 

(51) 

Of course, H„(—x) = ( — l) n H„(x), which implies that the exact 
determination of equation (51) depends on whether n is even or odd. 
It is straightforward to show that 

v(2n) = v(2n - 1) 

+ (-1)" sin (6) -*. exp (-^^T*.™^ (52) 
and 
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v(2n + 1) = v(2n) - (-1)" cos (6) -^= exp {-f)[^) #»»-»(^) 

for n *z 1. 

It is known that 8 

# (x) = 1, H^x) * 2x 

and 

H n+1 (x) = 2ztf n (z) - 2nH n ^(x) (53) 

for n ^ 1, so that in principle, the sequence {C n } can be computed 
by using equations (47), (52) and (53). However, in practice the dif- 
ference-equations (52) and (53) are unstable and care must be taken 
in their implementation. We found it best to use the following approach: 
Let [x] denote the largest integer not exceeding x and let n! ! = 
n!/2 tn/21 [n/2]!. [So that (2n - 1) !! = 1-3 -5- ■ • (2n - 1).] Now, one 
can use equations (52) and (53) to obtain stable algorithms for 



1,(71) = 



TMn) , , ^ _ rv& 



and t 2 (n) = 



nil 



2 [n/21 

which are used to compute 

C n = h{n)U(n). 

3.3 Numerical Results 

To obtain insight into the effect of hard-limiting, we assume that X t 
is bandlimited with a rectangular power spectrum. That is to say 

Ri{t) = r 2ttWt f ° r I T ' < °° * (54) 

In this case, 9 

S.(j/W) = C S(j) + a^ £ CJJf/W), (55) 

for —oo < / < oo. The functions F„(\) are denned as follows: 9 

\ir, -1 < X < 
[0, otherwise, 
and for n ^ 2, 



F 1 (X) = f' -lOKL (56) 
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F n (X) = . 



rnr 



E (-1)' 



X I + n 



- k 



k\(n - k)\ 



g | X | < n, (57) 



[0, otherwise, 

where M(n, A) = [n + |A|/2], We use [a;] to denote the integer part 
of x. 

As M(n, A) increases, it becomes increasingly difficult to accurately 
sum the alternating series displayed in equation (57). Since 



' — («M-¥ 



(58) 



(see Ref. 9), we used the asymptotic approximation (58) whenever 
F„(X) was required for n ^ 15. The threshold 15 was determined 
empirically. We attempted to keep our numerical estimates of the 
various spectra accurate to one percent relative error. 

We considered three cases: low modulation-index (r = 0.1), mod- 
erately high index (r = 1) and high index (r = 5). The results are 
displayed in Figs. 2 through 4. In each case, we have computed the 
effect of changes in the relative clipping level a = b/T. 

3.3.1 Low Modulation-Index 

When r « 1, it is evident from Fig. 2 and equation (55) that the 
principal part of the spectrum is well approximated by 



S.(j/W) 



Co *(/) + ~^ Ftf/W), for \f\<W, (5Q) 
0, otherwise. 



Using results from the preceding section, we can show that when 
r « 1, C, « r 2 exp (-T 2 ) (erf (a/yJ2)) 2 . Hence, for r « 1, 

s.(f/w) 

» k M + (erf (^)) 2 f2 6X 2 P ^- r2) FMW), \f\<W, (flQ) 

0, otherwise. 

That is to say, in the low-index case (T <<C 1), hard-limiting causes 
the principal part of the spectrum to be scaled down by the multiplica- 
tive factor (erf (a/V2)) 2 . 

Figure 2 illustrates the actual behavior [as opposed to the approxi- 
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2 4 6 8 

f/W 

Fig. 2— Spectrum of low-index PM wave (r = 0.1) for various relative clipping 
levels, a is the ratio of the clipping level b to the rms modulation-index r. 

mation (60)] of S v (f/W) for f > W. Notice that the tails of the spec- 
trum are raised considerably as the relative clipping level a decreases. 
This phenomenon is noticeable even for a = 4 (i.e., clipping at a 
four-sigma level). The increase in high-frequency content is appar- 
ently caused by the introduction of points at which the derivative of 
W(t) is discontinuous. 

3.3.2 Moderately High Modulation-Index 

The results for r = 1 are displayed in Fig. 3. The qualitative results 
are similar to those observed for r = 0.1 ; i.e., the principal part of the 
spectrum tends to decrease and the tails increase, as the relative 
clipping level a decreases. However, notice that for a particular value 
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of a, the frequency f/W at which the tails of the spectrum begin to 
rise is somewhat larger for r = 1 than for T = 0.1. 

3.3.3 High Modulation-Index 

When T = 5, the behavior of (V(t)) as a function of a is quite dif- 
ferent from that observed for T ^ 1 (see Fig. 1). In fact, small changes 
in a can change the discrete component of the spectrum, C = (V{t)) 2 , 
from a relative maximum to zero. Consequently, we originally suspected 
that the continuous part of the spectrum might change significantly 
as the discrete component changed from a relatively large value to 
zero. To check this point, four values of a were selected for testing; 




4 r / C 

f/w 



Fig. 3 — Spectrum of waves having moderately high modulation-index (r = 1) 
for various relative clipping-levels, a is the ratio of the clipping level b to the 
rms modulation-index r. 
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two at points of relative maxima for C and two where C was nearly 
zero. The results are displayed in Figs. 4a and b. 

From Figure 4a, we see that the qualitative behavior of S, (j/W) 
for high modulation-index (r = 5) differs significantly from that 
observed for T ^ 1. However, whether or not C vanishes does not 
seem to affect the general characteristics of the continuous part of the 
spectrum. Generally speaking as a decreases, the portion of the spectrum 
for < j/W < 1 tends to increase while the part for 1 < j/W < 10 
tends to decrease. Notice that a discontinuity is introduced at j/W = 1, 
a characteristic of lower modulation-index. Of course, limiting reduces 
the variance of the modulating signal X, , so that the discontinuity 
at j/W = 1 is not surprising. 

To examine the behavior of the tails of the spectra, look at Figure 4b. 
As in the other cases, decreasing a raises the tails of the spectrum, 
although the effect is considerably less than that observed for r ^ 1. 
In fact, for a > 3, we observed very little difference between the clipped 
and undipped cases. 

3.4 Comments 

Computation of the spectrum for large values of j/W, using the 
infinite series approach, is expensive since many terms are required. 
Moreover, accurate computation of C n and F„ , for large n, is difficult 
(if not impossible) because of numerical problems. Consequently, a 
good estimate of the spectrum for large j/W would be very desirable. 
Unfprtunately, we do not have such an approximation. However, it 
is possible to estimate the relative frequency j/W at which the tails of 
the spectrum for a clipped modulator will begin to depart from the 
undipped case. 

More precisely, consider equations (48) and (49) which constitute 
an algorithm for the computation of v{n). Since the tails of the spec- 
trum of W(t) can rise only when v(n) increases as a function of n, 
we need to examine v{n) for large n. 

For large n, we know that 6 

H 2n (x) = (-l) n 2 n (2n - 1)!! exp (x 2 /2)|cos (\An + 1 x) + oQjjj 

(61) 
and 

H 2n+i (x) = (-l) n 2 n+ *(2n - l)!!V2n+ 1 

• exp (z 2 /2)|sin ( V4n + 3 x) + O^j J \ • (62) 
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Hence, for large n we have the approximation 
v(2n) « v(2n - 1) - 2T sin (b)-yj- cos (aV2n - 3/2) 

-(- a r)n(^). m 

K2fl + 1) « v(2n) + cos (6) J- sin (aV2n + 3/2) 



f « \ V2n+ 1 ff ( 2k - l \ ,. .. 

•expl-r ^ II — F2— • (64) 



It is evident that the sequence /x n = IIlUi f( 2fc — !)) r ] I s a decreasing 
function of n for n < n = [r 2 -f- 3/2] and is an increasing function 
for n > n . For any finite a, m„ will ultimately exceed exp (a 2 /4) and 
the forcing function for v(2n) will begin to grow without bound. It 
follows that v(2n) will exhibit instability when n is so large that /z„ >>> 
exp (a 2 /4). This will certainly be true when (2n — l)/Y 2 ^ exp (a!/4) 
or for 

n^ |[r 2 exp (a 2 /4) + 1]. (65) 

An inspection of Figs. (2), (3) and (4) shows that equation (65) gives 
a good indication (at least for those cases considered) of the value 
of f/W where the tails of the clipped spectra begin a significant depar- 
ture from the undipped case. Equations (63) and (64) also indicate that 
for large n, v(n) will oscillate as it increases in magnitude. Consequently, 
it appears that it will be difficult to get tight bounds on S v (f/W) for 
large j/W. It is interesting that the carrier power is most significantly 
altered by the hard-limiting when r » 1 while the effect (of hard- 
limiting) on the tails of the spectrum decreases as V increases. 

IV. CONCLUSIONS 

In principle, equation (68) can be used to obtain either the auto- 
correlation function or the spectral density of any "reasonable" func- 
tion of a stationary gaussian process. We made use of the identity to 
study the spectral density of a sinusoidal wave which is phase- 
modulated by a hard-limited gaussian process. 

The main disadvantage of this approach (i.e., using an infinite series 
solution) is that numerical techniques are usually required to obtain 
a solution. Consequently, estimates of the tails of spectra (of interest 
for PM systems having high modulation-index) are difficult and 
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expensive to compute accurately. However, when no better alternative 
is available, our approach can be used to obtain numerical results in 
a relatively straightforward manner. 
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APPENDIX 

A Relation 

Let X t denote a stationary gaussian process with probability 
density function 

p ® = 2rT* GXP V 2T7 ' " °° < X < °° " ^ 

The joint probability density function for X (+T and X f is 

/ v 1 [ x\ - 2rx,x 2 + xf \ 

P(Xl ' Xi) = 2,rVl-r 2 GXP L" 2I-(l-r 2 ) J ' 

- oo < x, < * , (67) 



where 



R,(t) = (X l+r X,), T 2 = RM and r = && 



for — oo < r < oo . 

If G(x) is an exponentially bounded complex-valued function of the 
real variable x (i.e., there are real numbers u and v such that 16^(^)1 ^ 
v exp (u \x\) for — oo < x < oo) and if G(x) is of bounded variation 
on finite intervals, then we shall show that 



/OO /»00 

/ G(x l )G(x 2 )p(x l , x 2 ) dx t dx 2 



-t^F\L^^h 



(68) 



A.i Comments 

All the improper integrals displayed in this Appendix are defined 
in the sense of principal value. For example, equation (68) is more 
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precisely expressed as 



lim / / G(xi)G(x 2 )p(xi , x 2 ) dx x dx 2 



= ±mM\ lim f Gix) (£ Mx) ) dx 

„ = HI I 5T-.M J-T \0-X ' 



(69) 



Since G(x) is exponentially bounded and integrable over (— T, T) 
for all T ^ 0, it follows that the integrals in equation (69) exist and 
are bounded. Consequently, the improper integrals in equation (68) 
also exist. Moreover, one can show that 

lim / G(xi)G(x 2 )p(x 1 , x 2 ) dx y = f G{x^)G{x 2 )p{xi , x 2 ) dx y 

uniformly for — °o < x 2 < °° . 

It follows that the double limit in equation (69) can be expressed 
as an iterated limit and further that the improper integral in equation 
(68) can be evaluated as an iterated improper double integral. 

A.2 Proof of Identity 

We first establish equation (68) for all functions G(x) which are of 
bounded variation and satisfy 

f | G(x) | dx < oo . 

J-00 

In this case (following Rice, 2 Middleton 3 and Bennett 4 ), G(x) has a 
Fourier transform 

9(co) = f G(x) exp (-jo>x) dx (70) 

such that for almost all x, — oo < x < oo, 

G(x) = ^ /_" g(«) exp (jux) du. (71) 

Using equation (71), we can write 

/OO rtOO 

/ G(x 1 )G r (x 2 )p(a; 1 , x 2 ) cfo, dx 2 
•00 •»— 00 

= 47 i_ i_ \i_ ^^ exp ^* a ^ **7\i S ^ exp ^ 2 ^ *"*/ 

■p(^i , z 2 ) rfx! dz 2 . (72) 
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Since G(x) is of exponential order, one can show that all combinations 
of the integrals in the right-hand side of equation (72) converge uni- 
formly so that we are justified in interchanging the order of integra- 
tion to obtain 

/oo <•« 
/ G(x 1 )G(x 2 )p(x i , x 2 ) dxy dx 2 
•00 ■»— 00 

- i? L L swsw 

• ( / / p(x, , x 2 ) exp [j(oiiXi + co 2 z 2 )] dx x dx 2 \ dui du 2 . (73) 
Using the fact that 

/OO ilOO 

/ p(xi , x 2 ) exp [KuiXi + u 2 x 2 )] dxt dx 2 

= exp s — £■ (fa] ■+- 2ru l u) 2 + w 2 ) ) 

(see Ref. 5, pp. 30-35), and recalling that £n-o x*/n\ converges uni- 
formly to exp (x) for - oo < x < » , we obtain the following from 
equation (73): 

/CO f*00 

j GixJGixJpix! , x 2 ) dx Y dx 2 
-00 •'—00 

- t ^ fe L <*>"«" exp Hr -0 4 

• {^ f (ja, 2 ) n 9(o; 2 ) exp (-y co*) dco 2 j- (74) 

An application of Parseval's theorem yields 

£ /_" 9W(*i)" exp (-£ a, 2 ) rf„ = /_" <?(-*)(j£ P(x)) <** (75) 

and a similar expression for g and G. Using equation (75) in equation 
(74), we obtain equation (68). 

Now, we only assume that G(x) is exponentially bounded and of 
bounded variation on finite intervals. For A; = 1, 2, • • • , define 

q fa = faC*). for -k^x^k, 
[ 0, otherwise. 
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It follows that G k (x) is of bounded variation and that /" M | G k {x) \ dx < oo 
for k = 1, 2, • • ■ so that the above results imply that 



/ / G k (x l )G k (x 2 )p(x 1 , x 2 ) dxt dx 2 



(76) 



for k = 1, 2, • • • . 
We know that 



j£_ ( v 1 

dx n pw V2^r 5 V v^ W " n V\^ r 



HA-£-. ) exp (-^, 



[H„(z) is the Hermite polynomial of degree n], so that the series in 
equation (76) can be expressed as 



S k (r) = 2 a n , k r" 



where 



«„.* = 2^! 



Making use of the inequality 



7= [ G k (V2 Tx)H n (x) exp (-a: 2 ) dx 

/T •'-oo 



(77) 
(78) 



| # n (z) I < £ exp (f-)2 n/2 vnl, £ « 1.086435, 
and the assumption 



I G(x) | ^ v exp (m | a: |), — °o < x < oo, 
we can see that 

^ a n . t < (-§= ( exp (u\/2 Tx - x 2 ) dx) < oo , 

for all n and all fc. Hence, when | r \ < 1, Weierstrass' Tlf-test implies that 
the series S k (r) converges uniformly for A; = 1, 2, • • • . Consequently, 
we have 

/OO /»00 

/ G k (x l )G k (x 2 )p(x 1 , x 2 ) dxx dx 2 
•00 "—00 

for all t such that R x (t)/T 2 ^ ±1; i.e. for | r \ < 1. 



(79) 
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One can also show that 

/CO A 00 

/ G k (x y )G k {x 2 )'p{x l , x 2 ) dx x dx 2 
-00 J-00 

= f f G(x 1 )G(x 8 )p(x 1 , x 2 ) dx t dx 2 , (80) 
and that 

>™ £ °>®(£ p{x) ) dx = £ G{x)l S **) dx - (81) 

Substituting equations (80) and (81) into equation (79) yields the 
identity (68) for all r such that R x (t)/T 2 ^ ±1. 

We saw above that our identity can also be expressed as 



/ / G(x l )G(x 2 )p(x 1 , x 2 ) dXi dx 2 = 2_) a " r "> 

J-„ J-co »-o 



< 1, (82) 



where 

1 



-^= f Giy/2 Tx)H n (x) exp (-a: 2 ) dx 



(83) 



Since the left-hand side of equation (82) exists for r = ±1 and since 
a n ^ for all n, it follows that equation (82) is also valid for r = ±1. 
(The proof of this result is a Tauberian theorem. For example, see 
Ref. 10, page 427, exercise 13-34.) 
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